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The cross section for resonant ion-pair formation in the colhsion of low-energy electrons with 
HF"*" is calculated by the solution of the time-dependent Schrodinger equation with multiple coupled 
states using a wave packet method. A diabatization procedure is proposed to obtain the electronic 
couplings between quasidiabatic potentials of '^E''' symmetry for HF. By including these couplings 
between the neutral states, the cross section for ion-pair formation increases with about two orders 
of magnitude compared with the cross section for direct dissociation. Qualitative agreement with the 
measured cross section is obtained. The oscillations in the calculated cross section are analyzed. The 
cross section for ion-pair formation in electron recombination with DF^ is calculated to determine 
the effect of isotopic substitution. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

In low-tcniperaturc plasmas there is a significant con- 
centration of molecular ions. Processes such as disso- 
ciative recombination (DR) and resonant ion-pair (RIP) 
formation, that modify the charge and energy balance in 
these plasmas, are of great importance. To model the 
plasma environments correctly, a deeper understanding 
of the mechanisms and cross sections of these processes 
is needed. 

In DR, an electron is captured by the molecular ion 
losing its energy either to electronic or to ro-vibronic ex- 
citations of the resulting neutral molecule, corresponding 
to the "direct" or "indirecf'processes respectively of DR 
as proposed by Bardseley The molecule then stabi- 
lizes by dissociating into fragments. In the DR process, 
the fragments are neutral species, whereas in the RIP 
process the fragments consist of an ion-pair. 

The cross section of ion-pair formation in electron re- 
combination have been measured for a limited number of 
ions such as HD+ H+ %MM^ N0+ [|, 0H+ 

and IIF+ J,, Ji]J . Several of these measurements used the 
ion-storage ring technique 0, i, i, 0, S, i, 111. The ad- 
vantage of using storage rings is that the ion has time 
to relax into its vibrational ground state prior to the ex- 
periment, so there is a well defined initial state. Also, 
the ionic fragments can be separated from the neutral 
fragments in the bending magnets, and hence the final 
ion-pair channel is also well defined. Thus, theoretical 
studies on the ion-pair process can be compared to ex- 
perimental studies and give more insight into the physics 
of the reaction. 

So far only a few theoretical studies on resonant ion- 
pair formation have been presented. Only the systems 



HD+ I, [U, HeH+ [13, and H+ 0, [H, [11 have been 
examined. In these studies semi-classical methods [1, [11] 

or wave packets [3, [n m, m have been used to de- 
scribe the dynamics. The modeling of the RIP process 
is a theoretical challenge. Assuming, the direct process, 
a resonant state is created when the electron is captured 
by the molecular ion. Both the resonant state potentials 
and their corresponding autoionization widths must be 
well described. When the potential of the resonant state 
has crossed the ionic ground state potential, autoioniza- 
tion is no longer possible and the resonant state becomes 
electronically stable. It will then interact with the mani- 
fold of Rydberg states with potentials situated below the 
ion. To be able to capture the true dynamics of the RIP 
process, accurate potential energy curves, autoionization 
widths and electronic couplings need to be calculated 
over a large range of internuclear distances. The nuclear 
dynamics have to be explored from the Franck-Condon 
region where the electron is captured and into the asymp- 
totic region. The system may take different pathways on 
the way to dissociation into the ion-pair channel. Sev- 
eral interesting quantum effects may be studied, such as 
interferences due to the separate routes to the ion-pair 
limit yielding (Stiickelberg) oscillations J^] in the cross 
section for the reaction or tunneling through barriers in 
the potentials yielding resonant structures (shape reso- 
nances) in the cross section. 

In this paper we present calculations on the ion-pair 
formation in electron recombination with IIF+, i.e., the 
process: 



HF"" 



(1) 
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The cross section for this reaction has been mea- 
sured [9| for collision energies ranging from 0.0001 to 
1 eV, using the ion-storage ring CRYRING, where the 
F~ fragments were detected. Since the electron affinity 
of F is large, this ion-pair limit has a threshold energy 
of only 0.017 eV relative to the w = level of the ion. 
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This is within the rotational energy spread of the HF"*" 
target ion and therefore the measured cross section did 
not reveal any threshold effects. The ion-pair cross sec- 
tion was found to be relatively large, about 14 % of the 
DR cross section at a collision energy of 0.02 eV. The 
cross section for RIP shows interesting structures. Using 
photo ion-pair experiments, where the dissociation into 
the ion-pair -I- F~ is studied using photon excita- 
tion of the ground state X^Y,^ HF molecule into one of 
the bound Rydberg states, similar structures have been 
observed [H, [13, El] ■ Very recently, this RIP process 
has also been studied experimentally using the TSR ion- 
storage ring In this experiment, the absolute cross 
section was not measured. However, the shape and the 
structures in the relative cross section were very similar 
to those measured in CRYRING. Using an imaging tech- 
nique where the neutral H(n = 2)-t-F(^P3/2) fragments 
were detected, the rotational distribution as well as the 
population of the llF+{X^U3/2) and HF+(X2Hi/2) com- 
ponents were explored [l0| . 

We have previously studied the direct mechanism of 
DR of HF+ [lH. We calculated 30 resonant states of 
HF and a time-independent method was used to deter- 
mine the total cross section for dissociative recombina- 
tion. Autoionization from the resonant states was in- 
cluded in the model. The electronic couplings between 
the neutral states were neglected. It was thus assumed 
that the flux captured into a resonant state would follow 
that state diabatically out into the asymptotic region, 
without a redistribution of the flux. This calculation 
produced a total cross section for DR of HF+ in qualita- 
tive agreement with the measured cross section. Sharp 
threshold effects were seen where the asymptotic states 
become energetically open. The measured cross section 
below 0.04 eV was not reproduced and this was believed 
to be due to the neglect of the electronic couplings be- 
tween the neutral states. The lowest resonant state, 
in the quasidiabatic representation, dissociates into the 
ion-pair limit, H"*" -f F~. The calculated cross section 
for the direct dissociation into the ion-pair state was two 
orders of magnitude lower than the experimental results 
from CRYRING [§]. In the previous study, we presented 
the hypothesis that inclusion of the electronic couplings 
in the model might increase the cross section of ion-pair 
formation. Flux captured by higher resonant states could 
couple into bound Rydberg states that then predissociate 
by the electronic couplings with the ion-pair state. 

In the present study, we propose a diabatization proce- 
dure that allows us to obtain not only the quasidiabatic 
potentials, but also the electronic couplings among the 
neutral states. This method is applied to diabatize the 
^T,'^ states of HF relevant in dissociative recombination 
and especially ion-pair formation. The cross section is 
then studied by propagating wave packets on coupled 
states. Section [H] describes the calculation of relevant 
potentials and autoionization widths. The diabatization 
procedure is outlined in section [1111 while section |IV] de- 
scribes the treatment of the reaction dynamics. The re- 



sults are presented in section |Vl Atomic units are used 
throughout the paper unless otherwise stated. 



II. POTENTIAL ENERGY CURVES AND 
AUTOIONIZATION WIDTHS 

In our earlier study on HF potential energy 

curves and autoionization widths were calculated ab ini- 
tio by combining multi reference configuration interac- 
tion (MRCI) structure calculations with electron scat- 
tering calculations. In the present study, we are only 
interested in the ^T,~^ symmetry containing the ion-pair 
state. The potentials and autoionization widths of the 
resonant states situated above the ionic ground state po- 
tential are taken from the earlier scattering calculations. 
However, electronically bound states are recalculated us- 
ing the MRCI method on a much finer grid for internu- 
clear distances ranging from 1.0 ao to 9.0 ag. This is 
done in order to accurately resolve the avoided crossings 
among the adiabatic states induced by the couplings be- 
tween resonant and Rydberg states. For computational 
details we refer the reader to the previous study [l^ . We 
only present a brief summary of the calculation method. 
For the MRCI, the molecular orbitals are first determined 
using a SCF calculation followed by a MRCI calculation 
on the neutral ground state to calculate natural orbitals 
and hence a more compact representation of the orbitals 
is obtained. These natural orbitals are then used in order 
to carry out the MRCI calculations on ionic ground state 
(X^H) and the excited states of HF of symmetry. 
The MRCI calculation has reference configurations cre- 
ated by all excitations among the seven natural orbitals 
with highest occupations (3 cr-orbitals and 4 7r-orbitals) , 
except for the lowest (Icr) core orbital that is kept dou- 
bly occupied in all reference configurations. Single exci- 
tations from this set of reference configurations into the 
virtual orbitals are included. At each internuclear dis- 
tance, 25 roots are calculated. 

The same target wave function is used in the elec- 
tron scattering calculations carried out with the complex 
Kohn variational method [13] ■ By fitting the eigenphase 
sum of the transition matrix to a Breit-Wigner form [2lj 
both the autoionization with T(R) and the resonance en- 
ergy Eres{R) can be determined. The potential energy 
curve of the resonant state is obtained by adding the 
ionic potential to the resonance energy. In the present 
calculation, three resonant states of symmetry are 
included. 



III. DIABATIZATION PROCEDURE 

Rydberg states have the same character as the ground 
state of the ion [(lCT)^(2CT)^{3CT)^(l7r)^] plus an outer 
electron in a diffuse orbital. The resonant states are Ry- 
dberg states that converge to electronically excited ionic 
cores. These cores all have the (3f7) orbital singly ex- 
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cited. By following the configurations of the resonant 
states we obtain an "initial guess" of the quasidiabatic 
potential curves. This was the method applied to ob- 
tain the uncoupled diabatic potentials in our previous 
study ^9i- As mentioned above, the cross section for 
ion-pair formation, calculated by only including the res- 
onant state that diabatically correlates with the ion-pair 
limit, was two orders lower in magnitude than what has 
been measured Q. In order to improve this calculation, 
the electronic couplings between the neutral states must 
be included. Here, we propose a method where the adia- 
batic potentials and the "initial guess" of the quasi dia- 
batic potentials are used in order to determine not only 
the quasidiabatic potentials but also the electronic cou- 
plings. 



A. Extended two-by-two transformation 

Assume two adiabatic states interact in the vicinity of 
an avoided crossing. In the two-state model, the adi- 
abatic potential matrix V with matrix elements Vij = 
Vi{R)5ij can be a transformed into a diabatic potential 
matrix U using the transformation matrix P. The diag- 
onal elements of U are the diabatic potentials and the 
off-diagonal elements correspond to the electronic cou- 
pling. The transformation is given by 



U = P^VP, 



(2) 



where the unitary, orthogonal transformation matrix, in 
the two-state problem can be written as 



, cos[7(i?)] -sin[7(i?)] \ 
' sin[7(i?)] cos[7(i?)] i 



(3) 



Asymptotically, far from the avoided crossing, we assume 
that the adiabatic and diabatic potential curves are iden- 
tical. This assumption forces the rotational angle 7(i?) 
to be a function that goes from to 1. 

The two state model is now extended to include more 
states by treating each crossing as a two-by-two transfor- 
mation Fi , a rotation where only two states are included 
and the rest are left unchanged. If there are n crossings, 
the diabatic potential matrix, U, is given by 



U 



P^ip-iVPiPa.P^ 



(4) 



where the matrices P^ are of the form (in the case of 
rotation among state 1 and 2) 



Pi 



/ cos(70 
sin(7i) 





-sin(7,) 

cos(7i) 

1 

1 



(5) 



The rotational angles are assumed to have the following 
analytical form 



7,(i?) = -[l-Ktanh(a,(ii'-i?0)]. 



(6) 



An analytical expression for the total transformation ma- 
trix is now obtained and this matrix depends upon the 
unknown parameters of the rotational angles. The total 
transformation matrix can be shown to be an orthogonal 
matrix. After the diabatic potential matrix, U, is set 
up we perform an optimization procedure, where the un- 
known parameters ai and Ri of the rotational angles are 
determined in order to optimize the agreement between 
the diagonal elements of U and the estimated diabatic 
potentials obtained by tracking the configurations as de- 
scribed above. When the parameters are optimized, also 
the electronic couplings are obtained as the off-diagonal 
elements of the matrix U. By construction [see equation 
(HI)], the obtained quasidiabatic potentials have the adia- 
batic potentials as eigenvalues. However, we do not show 
that the couplings elements of the nuclear kinetic energy 
operator (non-adiabatic interactions) for the quasidia- 
batic states are identically zero. Therefore, the proposed 
method do not produce any strict diabatic states : 2^ , 
but rather quasidiabatic potentials and couplings. 

Figure [1] shows a schematic picture of this transforma- 
tion in the case of three states and two crossings. 



B. Quasidiabatic potential transformation 

The diabatization procedure described above is now 
used to obtain the potential energy curves of HF of 
symmetry. The experimental cross section for ion-pair 
formation is measured for collision energies up to 1 eV . 
Therefore, only the three lowest resonant states are in- 
cluded in the diabatization procedure. More energy is 
needed in order to capture the electron and form the 
higher resonant states. Furthermore, from the structure 
calculations the ground state of HF and five Rydberg 
states are obtained. 



1. Model I: 8 coupled states 

In the first model, the ground state of HF 
labeled here with Un) is excluded in the diabatization 
procedure. The ground state is assumed not to couple to 
the excited states of HF. Hence, in model I eight states 
are included in the diabatization procedure, five quasidia- 
batic Rydberg states [labeled with U22, ■ ■ ■ U^q] and three 
resonant states [C/77, Uss and Ugg]. State Ujr corresponds 
to the resonant state associated with the ion-pair limit. 
The optimization procedure is performed for internuclear 
distances 1.0 ao < i? < 9.0 ap. 

In Figure m the adiabatic (dashed) and quasidiabatic 
(solid) potential curves of HF are shown as well as poten- 
tial (thick dashed) of the HF+ ion. The potential curves 
of the quasidiabatic resonant states are black online. 

The calculated electronic couplings between the neu- 
tral states of HF are displayed in Figure[3l The electronic 
couplings between the Rydberg states and the resonant 
states [L/77 in (a), Uss in (b), and Ugg in (c)] are shown. 
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FIG. 1: (Color online) Schematic picture of the extended transformation. In (a), the three adiabatic potentials are shown. In 
(b) state 1 and 2 are rotated and finally in (c) a second rotation is performed among state 3 and 2. 
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FIG. 2; (Color online) Model I. Adiabatic (dashed) and qua- 
sidiabatic potential curves [solid (color online) lines] in the 
region of avoided crossings. Also the ground state of the HF^ 
ion is displayed with the thick dotted curve. The energy scale 
is relative to the ground vibrational state of the ion. 



Note that the electronic couplings are localized to the 
regions of the avoided crossings between the adiabatic 
states. 

The potential curves and autoionization widths are ex- 
trapolated to larger and smaller internuclear distances. 
The extrapolation toward smaller distances, is carried 
out as described in our previous study (Toj . The asymp- 
totic limits of the potential energy curves have to be de- 
termined. Here the experimental energies of the asymp- 
totic limits [1^ are used. Since the spin-orbit coupling 
in the system is neglected, the mean values of the energy 
levels with different J for the F atom and mean values 
for states with equal principal quantum number n for 
the H atom are used for the calculation of the asymp- 
totic limits. Wigner-Witmer rules [2^ are then applied 
to determine the number of states of symmetry that 



TABLE I: Asymptotic energy limits of the resonant states 
included in the present study. The energies are relative to the 
V = Q vibrational level of HF"'" . 

Atomic States Molecular States Energy (eV) 

1 H(n = 1) + F(2p) [/77 (Model II) -10.1552 

2 H+ + F^ Urr (Model I), Uu (Model II) 0.0170 

3 H(n = 2) + F(2p) (722, 1788, C/99 0.0437 

4 H(n = 4) + F(^P) U33, t/44 2.5934 

5 H(n = 1) + F{^P) U55, Uee 2.8993 



goes to each asymptotic limit. In table I, the asymp- 
totic limits for states included in the present study are 
listed. The potentials that are associated with neutral 
fragments are assumed to have reached their asymptotic 
limits at i? = 50.0 ag. The ion-pair state has an asymp- 
totic Coulomb form 

V{R) = Vf^nai - I - (7) 

where Vfinai is the asymptotic energy limit and a = 
17.581 ag is the polarizability of 25]. We assume 
that the ion-pair state has this form at internuclear dis- 
tances i? > 20 ao and interpolate between the calculated 
ion-pair state and the asymptotic form of the ion-pair us- 
ing spline interpolation. The electronic couplings will go 
smoothly to zero outside the region where the potentials 
are diabatized. 



2. Model II: 9 coupled states 

In the literature, there has been a discussion about a 
possible change of character of the ground state of HF 
when the internuclear distance is stretched [1^ [13] ■ At 
small distances, the ground state has an ion-pair char- 
acter, while as larger distances it goes covalently to the 
lowest limit H(n = 1)+F(^P). Hence, there is a large 
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FIG. 3: (Color online) Model I. Electronic couplings between 
Rydberg states and the resonant states (a) U77, (b) Uss, and 

(c) (799. 



avoided crossing between the ground state of HF and the 
ion-pair state. Similar avoided crossings are possesse d by 
the alkali halides such as LiF, LiCl and Nal [H, [H, ■ 
In model II, we examine the influence of this avoided 
crossing upon the ion-pair formation. Therefore, the low- 
est quasidiabatic potential (Un) is in model II associated 
with the ion-pair limit, while the lowest resonant state 
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FIG. 4: (Color online) Model II. Figure (a) shows the adi- 
abatic potentials (dashed lines) and the quasidiabatic poten- 
tials (solid lines) of HF in the region of avoided crossings. Also 
the ground state of HF"*" is displayed with the thick dashed 
curve. Figure (b) shows the electronic coupling between the 
ground state and ion-pair state. 



{U77) goes to the ground asymptotic fragments (see ta- 
ble[I]). Nine states are then diabatized as described above. 
In Figure H] (a) the resulting adiabatic and quasidia- 
batic potential curves are displayed. In (b) the electronic 
coupling between the ground state and ion-pair state is 
shown. The electronic couplings between the higher ex- 
cited states of Model II are similar to the corresponding 
couplings in Model I. As can be seen in Figured (b), the 
coupling element between the states with Un and U77 
potentials is an order of magnitude larger than the other 
coupling elements. The size of this coupling element is 
proportional to the energy difference between the adia- 
batic curves at the avoided crossing between these states. 
This large electronic coupling indicates that the nuclear 
dynamics most probably will go adiabatically in this re- 
gion. 
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IV. REACTION DYNAMICS WITH WAVE 
PACKETS 

The cross section for dissociating into the different 
channels, particularly into the ion-pair state, is calcu- 
lated by propagating wave packets on the coupled neutral 
states. The electron capture will initiate wave packets on 
the resonant states. For resonant state i (with i > 7 in 
model I and II mentioned above), we have [3l| 

^,(t = 0,R) = ^l?^Xv=o{R), (8) 

where Ti is the autoionization width for the resonant 
state and Xv=o is the v = vibrational wave function 
of the ionic ground state. 

It is here assumed that the electron is captured by 
the molecular ion in its lowest vibrational state. This is 
adequate when we compare with measurements using the 
ion-storage ring where the ion vibrationally relax prior to 
the experiment However, the ions in the storage ring 
will be rotationally excited. Furthermore, the spin-orbit 
coupling is excluded in the calculation of the relevant 
potentials. The fine-structure splitting of the 
and components of HF+ is about 30 meV |33 | 

and both of these components will exist in the storage 
ring experiments [lo| . 

The wave packets are initiated only on the resonant 
states. The non-adiabatic coupling between the ioniza- 
tion continuum and the Rydberg states is thus neglected 
and therefore, only the direct mechanism of DR is stud- 
ied. 

The nuclear dynamics is explored by solving the time- 
dependent Schrodinger equation 

* R) = _i.i^Vj>(t, R) + U{R)^{t, R) (9) 

numerically by using a Cranck-Nicholson propagator [s^ . 
In this study we have propagated wave packets on both 
coupled and uncoupled potentials. In the uncoupled case 
U is a diagonal matrix containing only the quasidiabatic 
potentials of the resonant states. Autoionization from 
the resonant states are included by letting these poten- 
tials be complex when their potential curves are situ- 
ated above the ionic potential curve. In the local ap- 
proximation for treating autoionization, also called the 
"Boomerang model" [31|, the resonant state potentials 
are given by 

UuiR) = Uu{R)-i^^-^. (10) 

In our previous study of HF [l9| we examined the validity 
of this approximation, and indeed it is valid for the states 
included in the present study. 

When the wave packets have been propagated out into 
the asymptotic region they are projected onto energy nor- 
malized eigenstates of the fragments [$^(i?)] and the 



cross section can be calculated as 

a,(i?) = ^.g|(<I)^(i^)|vI/,(^^,i^))|^ (n) 

where g is the multiplicity ratio of the neutral state and 
the ionization continuum. In the present study a grid 
ranging from R = 0.8 ag to i? = 300 &o with grid steps of 
dR = 0.01 ao is used. The wave packets are propagated 
with a time step of dt = 0.1 a.u.. For the uncoupled po- 
tentials the final propagation time is t final = 1000 a.u., 
while t final — 4000 a.u. is needed for the coupled sys- 
tems. In the coupled state calculations, part of the wave 
packets are trapped in bound Rydberg states before they 
predissociate and therefore a longer propagation time is 
needed to reach convergence. 



V. RESULTS AND DISCUSSION 

The wave packets are now propagated on the coupled 
potentials of symmetry where three resonant states 
are included. As mentioned above, model I includes five 
Rydberg like potentials and model II also includes the 
ground state of HF and the electronic coupling between 
the ground state and the ion-pair state. To test the role 
of the electronic couplings, we also propagate the wave 
packets on uncoupled potentials. Furthermore, the roles 
of rotation as well as the spin-orbit splitting of the target 
ion state are addressed. Finally, the isotope effect on ion- 
pair formation is studied by calculating the cross section 
for the DF+ ion. 



A. Coupled system Model I 

When the wave packets are propagated on the cou- 
pled potentials, parts of the wave packets that are ini- 
tially captured into the resonant states will couple into 
the bound Rydberg states and temporary get trapped 
in these states before predissociation is induced by the 
electronic coupling between the bound state and states 
that are open for dissociation. Feshbach resonances are 
formed that creates oscillatory structures in the cross 
section for ion-pair formation. In order to obtain con- 
vergence of these structures, the propagation time must 
be long enough for the wave packets to form the tem- 
porary trapped resonant state and subsequently predis- 
sociate. We have found however that after a time of 
t = 3200 a.u. the structures in the cross section do not 
change significantly. Therefore at a final propagation 
time oit final = 4000 a.u. the calculation is converged. In 
order to avoid problems with reflection against the end 
of the grid or against the absorbing potentials for some 
of the states a much larger grid than in the previous 
calculation [l^ was used. Figure \E\ shows the ion-pair 
cross section calculated for different propagation times, 
where all 8 coupled states in model I are included. For 
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FIG. 5: (Color online) Time dependence of the ion-pair for- 
mation cross section calculated using the 8 couples states of 
model I. At t=4000 a.u, the cross section has converged. 



a propagation time of t — 200 a.u. only the direct disso- 
ciation along the ion-pair state is seen. This cross sec- 
tion underestimates the measurements with about two 
orders of magnitude. At t — 600 a.u. on the other hand, 
flux from the higher resonant states has predissociated 
through the Rydberg states into the ion-pair state caus- 
ing an increase in the cross section. However, the calcu- 
lation at t = 600 a.u. is not converged with respect to 
time. We can see that while a time oi t — 1000 a.u. is 
enough to obtain the amplitude of the cross section, at 
later times, structures appear and as mentioned above, 
these structures do not change significantly after a prop- 
agation time of t = 3200 a.u.. Longer time is needed in 
order to converge the structures at lower energies. The 
final propagation time is set to t final = 4000 a.u.. Such 
behaviour has been observed previously ■ 

In order to understand the oscillatory structures found 
in the ion-pair cross section shown in Figure [51 a test cal- 
culation is carried out where only the ion-pair state (C/77) 
and the third Rydberg state (L/44) as well as the couplings 
between these states are included. These potentials are 
displayed in Figure [H] (a) where also the ion potential 
(thick dashed) curve is shown. The dotted lines are the 
vibrational energy levels of the Rydberg state. The cal- 
culated cross section for ion-pair formation [see Figure [6] 
(b)] shows sharp dips and the positions of these dips 
match perfectly with the energy eigenvalues of the vibra- 
tional levels of the Rydberg state, marked with crosses 
in the figure. Note that this kind of structure disap- 
pears at the energies larger than the dissociation limit of 
the Rydberg state. For this test calculation where only 
two coupled states are included, the alternative time- 
independent approach used in our previous study 
is also applied in order to check the convergence of the 
wave packet propagation. A very good agreement for the 
cross section calculated using the two different methods 




Interaction energy (eV) 



FIG. 6: (Color online) In (a), the potentials included in the 
test calculation are displayed. These are the ion-pair state U77 
(red online) and the bound Rydberg state U44 (blue online). 
The energies of the vibrational levels of the Rydberg state are 
displayed with the dotted (blue online) lines. The potential 
of the ion is the thick dashed curve. In (b) the cross section 
for ion-pair formation obtained using the test calculation is 
shown. The structures in the cross section can be explained 
by temporary formation of the bound vibrational levels of 
the Rydberg state. The energies of these vibrational levels 
are marked with crosses. 



is found. 

In model I, five Rydberg states are included and each 
will contribute with these kinds of oscillations in more or 
less overlapping energy intervals causing the rich struc- 
ture seen in the ion-pair cross section for the coupled sys- 
tem. In order to study the effects on the ion-pair cross 
section from the couplings between the neutral states, 
one state at the time is added to the model. The re- 
sulting ion-pair cross sections are shown in Figure [T] 
The first calculation only contains the ion-pair state and 
as seen before, the cross section is underestimated by 
about two orders of magnitude. Then, the bound Ryd- 
berg states {U22, ■ ■ ■ Uee) are added one after the other 
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FIG. 7; (Color online) Cross section for the ion-pair state 
when different number of states are included in the coupled 
system 



and this causes the dips in the cross section mentioned 
above. Note that the magnitude of the cross section 
is not dramatically affected. When the first Rydberg 
states is added, the cross section decreases as expected 
and some structures emerge. Adding the second Ryd- 
berg state only affects the cross section slightly. This is 
expected since the magnitude of the electronic coupling 
between the ion-pair state and the second Rydberg state 
is about one order of magnitude smaller than the elec- 
tronic coupling between the ion-pair state and the first 
Rydberg state. When the third Rydberg state is added, 
more pronounced structures emerge. Adding the fourth 
and fifth Rydberg state does not effect the cross section 
significantly since the couplings to these states are rel- 
ative small. When the second resonant state (t/gs) is 
added to the model, the cross section of the ion-pair is 
significantly increased since flux is predissociating from 
this state through the Rydberg states and into the ion- 
pair state. This can be seen by the black dashed curve 
for the cross section in figure [71 It should be noted that 
at low collision energies, the second resonant states has 
the largest electron capture probability at low collision 
energies. Also the sizable coupling elements between the 
second resonant states and the Rydberg states explains 
the increase of the ion-pair cross section when this state 
is added. When the third resonant state is added, the 
cross section do not change significantly, except for the 
energy interval between 0.3 — 0.4 eV. This is the energy 
window where the crossings between this state and the 
Rydberg states occur, see Figure [2] 

In the present study, the indirect mechanism of disso- 
ciative recombination is neglected. We thus neglect the 
non-adiabatic coupling between the Rydberg state and 
the ionization continuum. Induced by this coupling, the 
electron could be directly captured into a ro-vibronically 
excited Rydberg state that then predissociate by the elec- 



tronic coupling to the ion-pair state. This could in princi- 
ple affect the cross section for ion-pair formation. How- 
ever, it is shown that the indirect mechanism is most 
important at very low collision energies, typically below 
0.1 eV [3^, [13, nil . We believe the inclusion of the indi- 
rect mechanism could affect the cross section at very low 
collision energies and furthermore it would cause sharp 
structures in the cross section. These structures are not 
resolved in the measured cross section using the stor- 
age ring experiment 0. The indirect mechanism is best 
adressed using Multi-Channel Quantum Defect Theory 
(MQDT) [lil. However, standard MQDT calculations 
do not give final state distributions and hence could not 
obtain the ion-pair cross section. The indirect mech- 
naism has been included in wave packet studies 'ss'l . This 
yielded accurate branching ratios, but the magnitudes of 
the partial cross sections were difficult to converge. 

The partial cross sections for all states included in 
model I are shown in Figure \8\ In (a) , the cross sec- 
tions for the resonant states are displayed. In (b) the 
cross section for first Rydberg state is shown, while in (c) 
also the cross sections for the higher Rydberg states are 
displayed. Only the lowest Rydberg states and the two 
resonant states are energetically open at lower energies 
(see Table I). The lowest Rydberg state has a barrier in 
its potential which explains the sharp onset around 0.4 
eV. The resonances below (and above) the barrier are 
so-called tunneling shape resonances. Similar shape res- 
onances have been found in for example the cross section 
for dissociative recombination of HcH+ [s^ . 



B. Model I vs. Model II 

In model II, the crossing of the ground state and ion- 
pair state are included and we investigate the importance 
of such a crossing on the dynamics. The cross section 
for ion-pair formation, calculated using the two different 
models is displayed in Figured Note that the inclusion 
of the curve crossing does not affect the magnitude of 
the cross section. The large electronic coupling between 
the ground state and ion-pair state (see Figure |4|) forces 
the wave packet to propagate adiabatically through the 
region of the crossing. At the higher energies, model II 
gives a smoother cross section. In model I, part of the 
wave packet will oscillate back and forth around the min- 
imum ion-pair potential and hence cause interferences 
that require a longer propagation time for convergence. 
In model II however, the possibility of dissociating into 
the ground state limit also exists, and less flux will os- 
cillate and interfere. This results in a smoother cross 
section. At the very low energies, the lack of total con- 
vergence with respect to time is more noticeable, and we 
believe that this is the main reason for the discrepancy in 
the cross section between the two models in that region. 
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FIG. 8; (Color online) Partial cross section for the 8 states 
included in model I, in (a) the cross sections for the resonant 
states are displayed, while (b) and (c) shows the cross section 
for dissociation along the Rydberg states. 



FIG. 9: (Color online) Cross section for ion-pair formation 
calculated using model I and model II of potential curves and 
couplings. 



C. Effects from spin orbit splitting and rotation 

In the structure calculation, the spin orbit coupling is 
neglected and therefore we do not see the splitting of the 
ground state for the ion into its two components 
and The splitting of these components las of 

the order of 30 meV [lH. In the experiments [1, ITol |. 
both spin-orbit components of the ion will be populated 
with some unknown distribution. The most recent ex- 
periment carried out using the TSR storage ring [lOj, 
shows fragments produced in dissociative recombination 
with the HF+ ion in the X'^I{y,f2 state with the rotation- 
ally excited J — 7/2, 9/2, . . . , 15/2 levels in the limit of 
vanishing collision energy. For the component, 
excitation of the rotational levels of J = 1/2, 3/2, . . . 13/2 
are found. These excited fragments of the ion will shift 
the threshold energy for ion-pair formation. The ion-pair 
limit has a threshold energy of 17 meV for the lowest 
Q. — 2>/2 state and J = 3/2 level. Note that even for 
rotationally relaxed ions, for the = 1/2 component the 
ion-pair limit will be energetically open. 

In order to simulate the effect of ion-pair formation in 
electron recombination with the X'^Ili/2 component of 
HF"*" , a calculation is carried out where the ion-potential 
is shifted upwards 30 meV relative to the neutral states. 
In figure [TU] the calculated cross sections are compared 
with the measured cross section from CRYRING [9!]. 
Since the ion potential is shifted upwards, less interac- 
tion energy is needed to reach the resonant state poten- 
tials and structures will therefore be shifted toward lower 
energies as can bee seen in the figure. 

Rotational excitation of the ion will also give rise to 
such effects. In order to calculate the ion-pair cross sec- 
tion of a rotationally excited ion such as HF"*" in X'^Ii'i/2 
state with J = 7/2, not only the zero-point energy is 
shifted but also an effective term J(J -I- l)/(2/ii?^) is 
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FIG. 10: (Color online) Cross section for ion-pair formation 
in electron recombination with a vibrationally and rotation- 
ally relaxed HF ions in the X^n3/2 [solid ( green online)] or 
states [dashed (red online)]. The cross section for 
electron recombination with HF^ in X^Il3/2 with J = 7/2 is 
shown with the dashed-dotted (cyan online) curve. The tri- 
angles show the measured cross section from CRYRING j^]. 

added to the diagonal elements of the potential matrix 
used in the wave packet propagation. The resulting cross 
section is shown Figure [TUl 

In the measurements there is a distribution of the f2 
components as well as the rotational excitation of the ion. 
This will cause the structures in the measured cross sec- 
tion to be washed out. This might help explain why less 
pronounced structures are found in the measured cross 
section compared to the theory. 



D. RIP of DF+ 

By changing the reduced mass of the system when the 
vibrational wave function is calculated as well in the nu- 
clear dynamics calculation, resonant ion-pair formation 
in electron recombination with DF"*" is studied. Again, 
model I is applied with 8 coupled neutral states. The 
cross section for both fl components of the ion, 
and X^Ili/2, are calculated and displayed in Figure [TT] 
The amplitude of the cross section of ion-pair formation 
in electron recombination with DF"*" is similar to that of 
HF+ , but the oscillating structures are in this case some- 
what smaller. The cross section for ion-pair formation 
with DF"*" has not been measured. 



E. Total cross section in symmetry 

The cross section for propagating wave packets on cou- 
pled states of ^S"*" symmetry is calculated using the two 
models I and II. In order to study the effects of the cou- 
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FIG. 11: (Color online) Cross section for ion-pair formation in 
electron recombination with a vibrationally and rotationally 
relaxed DF^ ions in the states or X'^Il-i/2- 
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FIG. 12: (Color online) Contribution to the DR cross section 
from the three lowest resonant states of ^E^ symmetry, cal- 
culated using uncoupled states, and the coupled models I and 
II. Also the measured cross section of DR 9] is displayed. 



plings on the DR cross section, in FigurefT2l the cross sec- 
tion for the coupled systems is compared with the cross 
section calculated by propagating wave packets on the 
three lowest uncoupled resonant states of symme- 
try. We also compare with the measured cross section 
for DR la]- 

When the electronic couplings are included in the wave 
packet propagation, the total cross section for dissocia- 
tion along the states is reduced in magnitude. For 
low interaction energies, flux is lost into Rydberg states 
that are energetically closed for dissociation, making the 
magnitude lower in the energy range 0.044 — 0.5 eV. 

The cross section below the sharp threshold around 
0.044 eV is somewhat increased when the couplings are 
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included. However, it is still several orders of magnitude 
smaller than the measured cross section. This might be 
explained by the existence of the excited fl — 1/2 and ro- 
tationally excited ions in the experiments. Also, the indi- 
rect process, induced by the neglected non-adiabatic cou- 
plings between ionization continuum and Rydberg states 
might be important here. 

Above 0.4 eV, there is no longer a large difference in 
magnitude of the cross section calculated for the uncou- 
pled and coupled systems. The couplings induce more 
structure in the total cross section. At even higher en- 
ergies more resonant states not included in the present 
model become important. In our previous study [l^, 9 
resonant states of symmetry were included. 

Also resonant states of other symmetries not consid- 
ered in the present study, play an important role in DR 
of HF+. For the symmetry, we noticed [l^ that the 
potentials as well as the partial cross sections were very 
similar to those of symmetry. The only large differ- 
ence was that the lowest state of '^E+ symmetry is as- 
sociated with the ground state limit asymptotically. We 
can therefore expect a similar behavior if we introduce 
couplings between states of '^E+ symmetry. The lowest 
states of and '^n play a less important role in DR, 
since the partial cross sections were already lower than 
the experimental DR cross section in the whole energy 
range. 

VI. CONCLUSIONS 

We have calculated the electronic states of ^E"'" sym- 
metry relevant for ion-pair formation in electron recombi- 
nation of HF+ using a combination of scattering calcula- 
tions and MRCI structure calculations. We then propose 
a diabatization procedure, where not only the quasidi- 
abatic potentials are obtained, but also the electronic 
couplings between the neutral states. The diabatization 
procedure is based on the construction of an orthogonal 
transformation matrix written as a product of successive 
two- by- two transformations. 

The nuclear dynamics are described using wave packets 



that are propagated along coupled system including up 
to nine states. Autoionization is included using the local 
approximation. With our models we obtain very good 
agreement in magnitude with the experimental cross sec- 
tion measured in CRYRING [9]. Inclusion of the elec- 
tronic coupling increases the cross section by about two 
orders of magnitude. Our calculated cross section shows 
structures similar to those in the measured cross section, 
but the structures are more pronounced and somewhat 
shifted relative to the experimental ones. 

The present study do not include spin-orbit coupling 
effects. To address these effects, the ion-pair cross section 
for electron recombination with the component 
of the ion, which lies about 30 meV higher in energy, is 
calculated. A cross section without a sharp threshold and 
shifted structures is obtained. Also recombination with 
rotationally excited ions (J = 7/2 is calculated) which 
caused shifted structures in the cross section and pro- 
duced a larger value of the cross section at low collision 
energies. 

The ion-pair formation in electron recombination with 
DF"*" has also been calculated using the same potential 
energy curves, widths and couplings, but with a different 
reduced mass in the dynamical calculation. 

By adding up the partial cross sections of the states 
of symmetry, included in the coupled system we 
examine the role of the couplings to the total cross 
section for DR of HF"*". We see clear indications of a 
reduced cross section for DR when electronic couplings 
are included. Better agreements with measured cross 
sections is obtained. 
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